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Abstract: Major facilitator superfamily (MFS) is a large class of secondary active transporters 
widely expressed across all life kingdoms. Although a common 12-transmembrane helix-bundle 
architecture is found in most MFS crystal structures available, a common mechanism of energy 
coupling remains to be elucidated. Here, we discuss several models for energy-coupling in the 
transport process of the transporters, largely based on currently available structures and the 
results of their biochemical analyses. Special attention is paid to the interaction between protona- 
tion and the negative-inside membrane potential. Also, functional roles of the conserved sequence 
motifs are discussed in the context of the 3D structures. We anticipate that in the near future, a 
unified picture of the functions of MFS transporters will emerge from the insights gained from 
studies of the common architectures and conserved motifs. 
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MFS Transporters 

Major facilitator superfamily (MFS) is the largest 
class of secondary active transporters and is present 
in cells across all life kingdoms. ! So far, over 15,000 
genes have been identified as encoding MFS trans- 
porters.? For example, 25% of prokaryotic transport- 
ers belong to the MFS family,? whereas in the 
human genome, genes coding for at least 110 MFS 
proteins have been identified. One of the defining 
features of transporters is that they commonly use 
the electrochemical potential of one substance, for 
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example, from ATP hydrolysis (for primary active 
transport) or an ion gradient (for secondary active 
transport), to drive the transport of another sub- 
stance (i.e., the substrate). Many MFS transporters 
utilize proton-motive force (PMF) to drive the trans- 
port process.” In such cases, an MFS transporter 
often consumes one proton to transport one molecule 
of an electroneutral substrate. Thus, MFS transport- 
ers appear to be more energy-efficient in terms of 
the stoichiometric ratio of substrate to protons, com- 
pared to the other major class of active transporters, 
the ATP-Binding Cassette (ABC) transporters. The 
ABC transporters consume two ATP molecules 
(equivalent to about 6 protons) per transport cycle. 
On the basis of transport direction of the substrate 
relative to that of the driving substance, MFS trans- 
porters can be classified into symporters, antiport- 
ers, and, in case of a driving substance being absent, 
uniporters. Considering the countless types of mole- 
cules that a cell encounters during its lifespan, the 
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Figure 1. Representative crystal structures of MFS transporters. LacY/1PV7, GlpT/1PW4, FucP/307Q, and YajR/3WDO are 
shown in ribbon diagrams. Cavity helices are colored in blue, rocker helices in cyan, and support helices as well as the amphi- 


pathic helices a6-7 in golden. 


number of transporter genes is quite small. There- 
fore, many transporters have to accommodate 
numerous types of ligands, and such poly-specificity 
is likely to be a common property of MFS 
transporters. ”S 

This review discusses common energy-coupling 
mechanisms in MFS transporters, and uses PMF- 
driven electrogenic transporters as an example. 
Many parts of the discussion are, however, applica- 
ble to other types of secondary active transporters. 


Two Major Conformations of the MFS 
Transporter 

An alternating access mechanism for substrate trans- 
port was proposed by Jardetzky in 1966 for transport- 
ers.? According to this mechanism, a transporter 
differs from a channel molecule in that a transporter 
does not open simultaneously to both sides of the 
membrane, whereas a channel molecule may. Thus, a 
transporter often has two major conformations, 
namely the inward-facing conformation (Cm) and the 
outward-facing conformation (Cout), in which the sub- 
strate is accessible to the cytosol/cytoplasmic and 
periplasmic/extracellular space, respectively. Such a 
large conformational change not only drives the sub- 
strate transport, but may also permit the transporter 
to transduce signals.!? The key question that is asso- 
ciated with the alternating access model and remains 
unanswered to this day is as to how external driving- 
energy is coupled mechanistically to the physical con- 
formational change of the transporter. 

Analysis of crystal structures of MFS transport- 
ers, for example, LacY (PDB ID: 1PV7)," GlpT/ 
1PW4," FucP/307Q,! and YajR/3WDO," revealed 
that these proteins possess a 12-transmembrane 
(TM) helix core consisting of two pseudo symmetri- 
cal domains, namely the N-domain (TMs 1-6) and 
the C-domain (TMs 7-12), and indeed have two 
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major conformations, that is, the Cm and Cout states 
(Fig. 1). Between the two domains there is a central 
cavity serving as both the binding site and transport 
path for the substrate. An MFS transporter usually 
consists of four, pseudo-symmetrical, three-helix 
repeats.’ The first helix in each repeat (i.e., TMs 1, 
4, 7, and 10) contributes to the formation of the cav- 
ity, and is thus referred to as the cavity-helix. The 
second helix in each repeat (i.e., TMs 2, 5, 8, and 11) 
is usually long and curved, possessing a banana-like 
shape. Together, these four helices form the side 
walls of the central cavity, and are referred to as 
rocker-helices, to indicate the fact that they are 
directly involved in inter-domain conformational 
changes. The third group of helices (i.e., TMs 3, 6, 9, 
and 12) are located at both ends of the longest 
dimension of the TM core. They usually do not 
directly contribute to the formation of the cavity. 
However, they do contact the lipid bilayer, thus 
being involved in hydrophobic mismatch interactions 
with the surrounding membrane, and are referred 
to here as support-helices. On average, the TM 
region of support-helices are usually about one-turn 
shorter than that of the other two groups.'” Thus, 
during the major conformational change, support- 
helices are likely to move less than the other two 
groups of helices do, relative to the membrane. In 
order to specify the membrane-related position of a 
given residue, the membrane-embedded part (about 
20 residues or five helix turns) of each TM helix are 
divided into five regions, and are numbered as 1—5 
starting from the periplasmic side (not from the N- 
to C- termini of the TM helix) (Supporting Informa- 
tion Fig. S2 in Ref. 17). For example, TM7.3 depicts 
the central region of TM7. The pseudo symmetrical 
architecture is a major structural feature of MFS 
transporters and forms the bases of our understand- 
ing of energy-coupling mechanisms. 
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On the bases of available crystal structures of 
MFS transporters, a rocker-switch hypothesis of the 
conformational change between the two major states, 
Cout and Cm, was proposed, representing the mecha- 
nism behind the alternating access model suggested 
earlier.” According to this mechanism, each the N- 
and C-domains maintain their structures as approxi- 
mate rigid bodies and change their relative posi- 
tions,! allowing the central cavity to access the 
cytosol side of the membrane in the Cm state, and to 
the periplasmic side in the Cout state. Nevertheless, 
rigid-body rotation based conformational change is 
not the only probable mechanism. Other types of 
non-MFS transporters may have different ways of 
using the alternating access mechanism.’ The 
rocker-switch mechanism has been investigated for 
MFS transporters that were solubilized in detergent. 
For example, double electron-electron resonance 
(DEER) has been used to monitor changes of inter- 
domain distances on both the periplasmic and cytosol 
side of the TM core, for example, in the sugar-H* 
symporter LacY from Gram-negative Escherichia 
coli” and in the drug-H* antiporter LmrP from 
Gram-positive Lactococcus lactis.*! In these experi- 
ments, changes occurring at the inter-domain dis- 
tance on the two sides of the membrane are often 
correlated reciprocally. Further, results from studies 
on chemical modifications, such as Cys-scanning (e.g., 
in LacY?”3), dicyclohexyl-carbodiimide (DCCD) reac- 
tion, and PEGylation (e.g., in MdfA**), also support 
the notion of a rocker-switch mechanism. Therefore, 
the “rocker switch”-like conformational alternation is 
a well-established mechanism that links the pseudo 
symmetry of the MFS structure to the alternating 
access mechanism of general transporters. 


Energy Sources Driving Conformational 
Changes of MFS Transporters 
Applying the rocker-switch model dictates that an 
MFS transporter functions by cycling its conforma- 
tion between the two major states. From an enzymo- 
logical point of view, the Cout conformation can be 
considered as the ground state of a PMF-driven, elec- 
trogenic, MFS transporter; and the Cm conformation 
can be considered as the excited state (Fig. 2). Transi- 
tion from the ground state to the excited state is 
energized by protonation inside the central cavity, 
whereas the transition from the excited state back to 
the ground state is more or less automatic, being 
driven by conformational energy stored in the excited 
state. A key question associated with the rocker- 
switch mechanism is how the two common functions 
of an MFS transporter are coupled, namely the proto- 
nation (or binding of other ions) and conformational 
alternation between the Cout and Cy, states. 

It is commonly accepted that in order for the 
MFS transporter to transport a substrate sustain- 
ably against its own concentration gradient (i.e., 
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Figure 2. Schematic diagram of the transport cycle of MFS 
protein. 


Akts] = RTin(Slr/[Sl_)>0, where subscripts “L” 
and “R” stand for the loading and releasing states), 
external driving energy is essential. In addition, in 
the case of a charge-carrying substrate moving 
against the membrane potential, the electrochemical 
potential of the substrate (including an extra electro- 
static energy term) must be compensated by an 
opposing term from the external driving energy. The 
total external energy of a PMF-driven transporter is 
Auy + Aur+)+ Aug,” and should be equal to the 
change of Gibbs free energy (AG, i.e. the negative 
value of the heat released to the environment) of the 
transport cycle. This free-energy change is necessar- 
ily to have a negative value according to the Second 
Law of thermodynamics. Within the total energy, 
the part of |Awy +Apy+)| is considered to be the 
driving-energy, and Apis; is the result of the trans- 
port cycle. In other words, part of the external 
driving-energy is used to move the substrate against 
its own electrochemical potential, and the rest is 
converted to heat. Thus, energy efficiency of the 
transport cycle may be defined as Aprg/|Auy 
+Aup+j| (< 1). At the micro-scale, an MFS trans- 
porter functions similarly to Maxwell’s Demon (in 
Theory of Heat, 1871), who senses a specific type of 
molecules (substrates) and opens a gate for them 
allowing a movement of the substrate against its 
concentration gradient. The entropy of the substrate 
decreases during the transport process, at an 
expense of part of the driving energy being con- 
verted to heat. The total entropy increases to satisfy 
the Second Law of thermodynamics. 

Despite of this simple thermodynamic notion, 
the nature of the mechanism behind this energy- 
coupling process has to this day remained a mystery. 
Specifically, how does the binding of a 1-Da proton 
drive the large conformational change of a 50-kDa 


Energy Coupling Mechanisms of MFS Transporters 


A coupled coupled coupled coupled 
4 
Court AG? oy; Apy 
= | 
2 
oO Q 
S 
o ATT AGa(S) Cout | 
D ON AG,(H*) Weel 
_ AG, 
LL s... E AGp(S) 
$ = J ` a J 
Cin Cout 
States of electrogenic antiport 
B coupled coupled coupled 
p a = 7 SS 
+ 
> 
D Apy Q 
o oe -164/ (AG%o) Cout CoS 
be \ Sout 
2 p = = CouSH* 
È AG,(H°) ‘a 
WL 
z 
, = 
s. v v = 3 
Cin Cout 
States of symport 
C coupled coupled coupled coupled 
CoutH+ 
AG? o Apy 
> | 
D Aty C,,SH* Q 
2 AGa(H") K$ AG? so -Aly 
S n ar Cos | 
‘Out 
besss AG,(S) 
® TN Sou 
oc AG,(H*) L(H") 
L seme AGE TER 
a - ON 7 7 


States of electroneutral antiport 


Figure 3. Schematic diagram of the energy landscape of MFS transporters. A: Electrogenic antiporter. B: Symporter. C: Elec- 
troneutral antiporter. The free-energy landscape plot describes the thermodynamic relationship between different states, without 
attention to kinetic issues. The plot must satisfy the First and Second Laws of thermodynamics. Horizontal lines represent 
states. Tilted lines represent transitions between states. Transitions associated with the proton are indicated in blue, those with 
the substrate in red, and those with AY in green. Subscripts L, R, D, and E stand for energy terms associated with loading, 
releasing, differential binding, and elastic, respectively. “I>O” and “O>I” stand for the Ci,-to-Cout and Cour-to-C;, conforma- 
tional changes, respectively. In principle, since the transport process cycles, choice of the starting point is arbitrary. In this 
sense the starting and ending states are identical, only being differed by the release of heat (Q) during one transport cycle. 
Thus, the end state must be below the starting state. Neighboring states may be coupled tightly. In such a case, their sequen- 
tial order may be arbitrary. Locally, any transition of positive AG must be driven by a neighboring transition of a negative AG. 


(Also see Appendices). 


transporter such as MFS. Figure 3 shows energy 
diagrams for both a symporter and an antiporter 
(also see Appendices A and B for detailed discussion 
on energy terms). For any step of an isothermal- 
isobaric transport cycle, a positive term of AG must 
be coupled to a negative term(s), to keep the sum of 
AG negative. 
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The AG terms associated with chemical concen- 
trations deserve some explanatory notes. The free 
energy change of a substrate, Ajs), can be divided 
into three terms: AG,(S), associated with substrate 
loading; AGp(S), what we call differential-binding 
energy between the loading and releasing states; 
and AGR(S), associated with substrate releasing. 
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AG,(S) determines how easily the substrate may 
bind to the transporter (e.g., from the periplasm for 
a PMF-driven symporter); and AGR(S) determines 
how easily the substrate may be released from the 
transporter (e.g., into the cytosol for a PMF-driven 
symporter). Both AG (S) and AGR(S) vary with the 
concentrations of the substrate on the two sides of 
the membrane, and may be coupled with energy 
terms of protonation/deprotonation. In contrast, 
AG»p(S) does not change with the substrate concen- 
trations on either side of the membrane, and thus 
is an intrinsic property of the transporter. In 
principle, AGp(S) is experimentally measurable (see 
Appendix C). A negative AGp(S) would indicate that 
the substrate-carrying conformational change from 
the loading state to the releasing state is thermody- 
namically favored, and thus it contributes in part to 
the driving force of the transport. The complex of 
the drug-H* antiporter MdfA with its substrate 
TPP* seems to fit this scenario? [see Fig. 3(C)]. 
There, substrate binding stabilized the substrate- 
releasing state more than the substrate-loading 
state. In contrast, if AGp(S) is positive, the binding 
affinity is higher in the loading state than in the 
releasing state. It may favor substrate loading and/ 
or releasing. Yet, in this case, external energy inputs 
are required to drive the conformational change, to 
overcome the positive AGp(S). The results of studies 
of the complex of the sugar-H* symporter LacY with 
its substrates seem to fit this scenario” [see Fig. 
3(B)]. There, substrate binding stabilized the 
substrate-loading state. Similar discussions can also 
be made for the proton associated free energy terms 
(i.e., Azz +y ete.). 

As a corollary of the above discussion, formation 
of an intra-molecular hydrogen-bond (H-bond) (or a 
salt-bridge bond) within the transporter may stabi- 
lize a particular state of the transport cycle, but can- 
not be a driving force for thermodynamic cycle of 
the transport. In fact, the same H-bond must break 
later in the process, and the energy released from 
the H-bond formation has to be paid back from other 
sources. Nevertheless, such short-term energy gains 
may influence the kinetics, for example the rates, of 
individual step(s). 


Protonation and Membrane Potential 


Membrane potential-driving hypothesis 

The plasma membranes of all living cells usually 
carry an electrostatic membrane potential (AY) of 
typically ~100-mV magnitude. If effective monovalent 
ions that generate this membrane potential were uni- 
formly distributed in a region of oth fraction of the 
membrane thickness on each side of the membrane, 
the required ion concentration would be estimated to 
be 5 + millimolar. For example, if the membrane 
thickness were 30 A and the ionic region were 3-A 
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thick (i.e., 6=0.1), the effective ion concentration 
would be 10 mM, a value that is biologically relevant. 
Such an ionic distribution is equivalent to 0.02 stand- 
ard charges (i.e., the charge of an electron or a pro- 
ton) per 1000 A? (an area value of the cross section of 
a typical membrane protein). Therefore, the mem- 
brane potential is a macroscopic physical property of 
the membrane. It is produced by fast Brownian 
motions of charge-carrying particles. Such clouds of 
charges on both sides of the membrane can be 
described best using an average-field hypothesis, 
which treats fast moving charges as a continuous dis- 
tribution, especially when considering the slow 
motions of the major conformational changes of a typ- 
ical transporter. 

Why is membrane potential important for PMF- 
driven transporters? From an energetic point of 
view, for a 100-mV membrane potential, the free 
energy of transferring a proton across the membrane 
(Aup= FAY, where F is the Faraday constant) is 
about —4 RT. In comparison, for a ApH of 0.6, corre- 
sponding free energy of transferring one proton 
(Aum+= —2.3 RTApH) is about —1.4 RT. Some 
MFS transporters that function under alkaline pH 
conditions may even perform transport under a neg- 
ative (i.e., acidic-inside) ApH.” In addition, pH 7.6 
(the typical pH inside a bacterial cell) is equivalent 
to ~15 effective free protons per 1-1m* volume, sug- 
gesting that a cell is necessarily subjected to fluctua- 
tion in pH. Since electrostatic membrane potential is 
universal to all living cells, it would have been bio- 
chemically and evolutionarily advantageous to select 
in favor of utilizing the membrane potential as part 
of the driving force for MFS transporters. In fact, 
the presence of a membrane potential increases the 
transport rate of LacY by a factor of over 20 (see 
Fig. 2 in Ref. 26). Because the relative dielectric con- 
stant (cr) of an aqueous solvent (~80) is much higher 
than that of a protein (~1), the electrostatic field 
(E = AV/(e,d)) of a membrane potential is focused on 
the protein part of the MFS transporter (with an 
effective thickness of d<30 A). In other words, in 
both Cout and Cy, states, the solvent-filled cavity of 
an MFS transporter does not bear a significant part 
of the electrostatic field. Because of this non-uniform 
distribution of the dielectric constants, an electro- 
static field jumps between the cytoplasmic half (in 
the Cout state) and the extracellular half (in the Cy, 
state) of the transporter during the major conforma- 
tional changes. Relative to the membrane, the posi- 
tional shift of the field may be even larger than the 
movement of charged particle(s) within the trans- 
porter. It is the movement of the charged particle(s) 
relative to the electrostatic field that determines the 
energy levels of a transporter. It is noteworthy, how- 
ever, that several reported crystal structures of MFS 
transporters have been presented in the so-called 
occluded conformation. !®?™?S Although such a 
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conformation may represent a transient intermedi- 
ate state of the conformational change, it is unlikely 
to be a stable state in the presence of membrane 
potential, and therefore should be treated as in vitro 
artifacts. If it were a stable state, the transport pro- 
cess would not be able to proceed. Taken together, 
Apy is a major part of the free energy input for a 
PMF-driven MFS transporter, especially for electro- 
genic transport. 

For PMF-driven transport, both a symporter 
and an antiporter get protonated in the Cout state 
and deprotonated in the C;, state. In the presence of 
a negative-inside membrane potential, a protonated 
site is necessarily subjected to an inward electro- 
static force.* This force is ~5 pN in strength (ie., 
~FAY/(é,d)), which is comparable with forces of typi- 
cal biological events at the molecular level.?°°1 
Here, we propose a simple, yet elegant hypothesis 
on how the major conformational changes in MFS 
are coupled to the free energy Amy (Fig. 4): Protona- 
tion provides a structural point for the electrostatic 
membrane potential to apply an inward force. 
Together with the force of a hydrophobic mismatch 
from the surrounding lipid bilayer, the electrostatic 
force generates a mechanic torque that promotes a 
relative rotation between the N- and C-domains, 
thus shifting the Cou-Cin equilibrium toward the 
Cy, state. The hydrophobic mismatch force that 
balances the electrostatic force is necessarily associ- 
ated with many small deformations on the protein- 
lipid interface (either membrane bending or inward 
shift of the transporter). On the one hand, the elec- 
trostatic force and the sum of hydrophobic mismatch 
forces are of equal size but opposite in their direc- 
tions, thus keeping the mass center of the TM core 
within the membrane (see Appendix D). On the 
other hand, the elongated shape of the MFS trans- 
porter (with the longest dimension defined by the 
distance between the two pairs of support-helices, 
that is, TMs 3, 6 and TMs 11, 12) ensures that, for 
each of the N- and C-domain, the two forces gener- 
ate a mechanical torque. Interestingly, the same 
architecture of a central cavity containing titratable 
residues inside, together with an elongated overall 
shape is also observed in some non-MFS transport- 
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ers, e.g., the MATE family of transporters.°??" In 


addition, some other non-MFS transporters form 
elongated dimers, as exemplified in the LeuT family 
member, Arg*/Agm antiporter AdiC,** and the 
NhaA family member, Na*/H* antiporter NapA.*° 
One explanation for such widely used 3D architec- 
tures is that the hydrophobic mismatch effect is 
essential for the conformational changes involved in 
transport function. Importantly, the membrane 
potential-driving hypothesis provides a unified foun- 
dation to explain common mechanisms of energy 
coupling with conformational changes of electrogenic 
MFS transporters, regardless of types of substrates. 

In a more general sense, the combination of 
electrostatic and hydrophobic-mismatch forces deter- 
mines the equilibrium position of a transporter (or 
any charge-carrying integral membrane protein) 
within the membrane.” Changing either the mem- 
brane potential or the charge distribution of the 
transporter results in positional changes of the 
transporter relative the membrane, which in turn 
may affect chemical properties of the membrane pro- 
tein and the kinetic equilibrium rates between its 
sub-states. For instance, in the MFS antiporter 
MdfA, a change in the electrochemical potential of 
H* shifts the equilibrium between the Cy, and Cout 
states in the presence chloramphenicol (an electro- 
neutral substrate).®” In fact, the interaction of a 
charge distribution with the membrane potential 
determines the orientation of an integral membrane 
protein in the membrane, and this is reflected in the 
widely accepted positive-inside rule (of charge distri- 
bution).*° ATP synthases and voltage-gated channels 
have been shown to utilize the membrane potential, 
either as a driving force or a means to regulate their 
gate functions.®**° Recently, membrane potential- 
dependent conformational change was proposed to 
be responsible for the sealing of the peptide- 
conducting channel in the SecYEG translocon, thus 
preventing ion leakage during translocation of sub- 
strate peptides across the membrane.*! Such a mem- 
brane potential-dependent, positional equilibrium of 
membrane proteins also implicates that, an in vitro 
study, be it functional or structural, would deviate 
from the in vivo situation to some extent, because of 
both the absence of a membrane potential and 
altered hydrophobic mismatch. 

Moreover, in an MFS transporter, not all Apy 
related energy is converted to heat in the Coy-to-Cim, 
step. Part of it may be stored in a form of elastic con- 
formational energy (AGr) in the transporter-membrane 
system, and is released during the Cy,-to-Cou 
conformational change. The stored elastic conforma- 
tional energy is the result of a balance between an 
external force (e.g., electrostatic force) and internal 
forces (e.g., those associated with a hydrophobic mis- 
match and other structural adjustments). In agreement 
with this elastic nature, the stored energy is roughly 
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proportional (albeit not necessarily linearly) to the 
external force. Withdrawing the external force (e.g., by 
deprotonation) would result in the release of the elastic 
conformational energy. Taken together, membrane 
potential-dependent positional and conformational 
changes of a PMF-driven transporter during the Cout- 
to-C;, conformational change are associated with elas- 
tic conformational energy stored in the Cy, state, which 
may serve as a driving force for the Cm-to-Cout confor- 
mational change. 


Protonation sites 

Proton-titratable residues are common features of 
PMF-driven MFS transporters, although their loca- 
tions may vary from one transporter to another. All 
reported crystal structures of PMF-driven MFS trans- 
porters contain such titratable residues inside their 
central cavities (Supporting Information Table S2 in 
Ref. 42), either acidic or histidine residues. For exam- 
ples, E. coli LacY, the most extensively studied pro- 
karyotic MFS transporter, contains a H322™1°?_ 
E3251TM10.3 pair inside the central cavity, which is con- 
served in many sugar porters and essential for the 
PMF-driven transport. In E. coli YajR (a putative 
antiporter), H225™7-3 and E320™™°3 reside close to 
each other in the 3D structure and are the only two 
titratable residues inside the cavity.'* E. coli YbgH (a 
dipeptide-H* antiporter) contains a conserved 
E21™! inside the central cavity, essential for proto- 
nation.*” While it is widely accepted that protonation 
is important for the transport process, the concept of 
interaction between protonation and membrane 
potential remains to be established. 

According to our membrane potential-driving 
hypothesis, it is not necessary for protonation to be 
fixed at one position during the Co,;-to-Cy, conforma- 
tional change, provided that there exists one of the 
titratable residues inside the central cavity assuming 
the protonation status at any given time. A proton 
may potentially be transferred between two (or more) 
binding sites during the conformational change, 
depending on difference of the pK, values of the 
titratable sites. For example, the multidrug- 
resistance antiporter E. coli MdfA contains conserved 
E26™!4 and D34™1 inside the central cavity, and 
both of them were shown to be required for transport 
activity.” On the one hand, in the Cout state, E26 is 
located at the apex of the cavity and buried more 
deeply than D34, suggesting a higher pK, for E26. 
Thus, E26 is likely to be the one that gets protonated 
in the Cout state. In the Cm state, on the other hand, 
D34 is located at the new apex of the central cavity 
and becomes buried more deeply than E26. Thus, D34 
becomes protonated in the Cm state in the absence of 
a substrate, while E26 becomes deprotonated. There- 
fore, during the Co,1-to-C;, transition, there appears 
to be a proton transfer occurring from E26 to D34, 
either directly (via a proton-relay path) or indirectly 
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(e.g., by a dropping-and-picking mechanism).’” Simi- 
larly, in the fucose-H* symporter FucP, a pair of acidic 
residues, D46™"? and E135™**, are located at the 
apex of the central cavity in the Cy, and Cout states, 
respectively, and flank the substrate from two sides.” 
Mutations of either D46 or E135 results in activity 
loss of fucose uptake. In addition, of the two acidic 
residues, E135 is directly involved in substrate bind- 
ing. Mutation of E135 abrogates the counter-flow 
activity (i.e., the ability of binding substrate), while 
that of D46 does not. These observations can be inter- 
preted as following: in the Cout state, substrate bind- 
ing from the periplasmic side induces the protonation 
of E135 by expelling solvent from the vicinity of E135. 
This proton is transferred to D46 during the Cout-to- 
Ci, transition. After releasing the substrate into the 
cytosol, D46 becomes solvent-exposed and thus depro- 
tonated. Therefore, proton transfer inside the central 
cavity during the Co,;-to-Cy, transition may not be an 
isolated phenomenon. 

For many non-MFS secondary active transport- 
ers, there exist two structural parts (or domains) 
that undergo conformational changes relative to 
each other during substrate transport.’ According 
to the membrane potential-driving hypothesis, the 
part that changes its electrostatic status in response 
to substrate binding is likely to be the one that 
moves more significantly, relative to the membrane. 
In agreement with such a prediction, the part that 
avoids charge-change is often involved in formation 
of an oligomer, which is likely to serve as the scaf- 
fold for the moving parts.?*3544 


Symporters, Antiporters, and Uniporters 

Given their common energy-coupling mechanism, 
how can PMF-driven symporters, antiporters, and 
uniporters be distinguished from each other? A 
major functional difference between a symporter and 
an antiporter is whether the Coy:-to-Cy;, conforma- 
tional change carries a substrate: Symporters do, 
but antiporters do not. Then, an interesting question 
arises as to whether it is possible for a given MFS 
transporter to function as a symporter for some sub- 
strates yet as an antiporter for others. If one proton 
was consumed for both uptake of one substrate and 
expelling of another per transport cycle, the process 
would be a PMF-driven exchanger, which is hypo- 
thetically possible but has so far rarely been 
reported. Moreover, if the influx and efflux were 
operating independently from each other, the trans- 
porter would allow a PMF-driven conformational 
change to occur in the absence of a substrate. Such 
a scenario would be equivalent to proton leakage 
and is likely to be detrimental for the cell. 
Thus, such a dual-functional MFS transporter is 
unlikely to exist. (Note, however, that E. coli KgtP 
transporter was reported to carry out both a- 
ketoglutarate uptake*® and arabinose efflux,*® albeit 
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a detailed mechanistic analysis is still missing.) Log- 
ically, a symporter must couple substrate binding 
with protonation, and an antiporter requires compe- 
tition of substrate binding with protonation (ie., 
coupling of substrate binding with deprotonation). 
These notions are supported by early experimental 
data for several MFS transporters”**’ and are likely 
to be generally true for most cation (H*/Na*)-driven 
transporters of different types.‘ In short, symport- 
ers and antiporters are fundamentally different from 
each other in terms of substrate-induced energy cou- 
pling. Therefore, it is important to understand the 
structural differences between symporters and 
antiporters. 


Symporters 
Since in symporters, protonation drives the Cout-to- 
Ci, conformational change, it is necessary that, for a 
symporter, protonation only occurs after substrate 
loading; otherwise, proton leakage might ensue. 
Under the condition that its substrate is present in 
limited amounts, a PMF-driven MFS symporter is 
likely to spend more time in its de-protonated Cout 
state, waiting for the binding of a substrate that 
moves from the extracellular/periplasmic side to the 
central cavity (unless an allosteric mechanism exists 
for other regulatory factors to trigger the Cy,-to-Cout 
conformational change). To couple substrate binding 
with protonation, an ideal symporter should have 
the following properties: The pK, values of the 
proton-titratable sites in Cou (ie., pKa out) in the 
absence of substrate binding is lower than pHout, 
preventing premature protonation. Substrate bind- 
ing from the outside increases pKa our to above pHout 
thus promoting protonation. In such a scenario, the 
free energy of substrate binding, AG ,(S)= 
—RTin((Slou/Kaut), may be partially used to raise 
the pKa out of the protonation site. In the Cm state, 
the pK, in value before substrate release is higher 
than pHj,, preventing premature deprotonation. 
Substrate release decreases pK,in to below pHi, 
thus promoting deprotonation. This paradigm 
implies that substrate loading in Cout is not later 
than protonation, and substrate release in Cy, is not 
later than deprotonation. In LacY, it is found that 
decreased affinity for sugar in the Cout state at alka- 
line pH is due specifically to an increased reverse 
rate,*® suggesting that, under normal pH, sugar 
loading promotes protonation which in turn stabil- 
izes sugar binding. In short, for a symporter, sub- 
strate binding is a prerequisite for proton binding. 
How is such a mechanism applied in a sym- 
porter structure? At least two types of implementa- 
tions have been observed in available structures of 
MFS symporters. One simple way is to change the 
micro-environment of the titratable residue by sub- 
strate binding, as exemplified in the above discus- 
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Figure 5. Coupling of substrate binding with protonation in 
symporters. Putative substrate binding site of E. coli YogH 
(PDB ID: 4Q65). A model substrate (alafosfalin) is included 
based on superposition of a homologous complex crystal 
structure (4IKZ). Selected residues that are potentially impor- 
tant for substrate binding are shown in stick models. A clus- 
ter of polar residues (Q18, E21, Y22, and K118) potentially 
plays a role of coupling substrate binding with protonation at 
E21. 


Furthermore, proton-titratable residues of symport- 
ers are often involved in a cluster of polar residues. 
Substrate binding may turn the cluster more electro- 
negative, either by attracting positive charges away 
from the cluster, and/or by pushing negative charges 
toward the cluster and thus, the cluster becoming 
more electronegative and easier to become proto- 
nated. For example, in the crystal structure of E. 
coli LacY/1PV7, residues Y236™7?, Dp240™™71 
R302TM9-1 K319™10-1 773997102 ang E325T™M10.3 
form a cluster of polar residues and are conserved in 
a number of sugar porters." This cluster may func- 
tion as a proton-relay network, in which the H322- 
E325 pair is the titratable site. In E325 mutant var- 
iants, both uphill (PMF-driven) transport and 
substrate-induced proton uptake are abolished, 
while substrate affinities may even be increased rel- 
ative to the wild-type form,®°°! suggesting that 
E325 is directly involved in protonation. In E. coli 
YbgH, a peptide-H* symporter, the titratable resi- 
due E21™1* is involved in a cluster of polar resi- 
dues, some of which directly interact with the 
negatively charged carboxy group of the substrate 
dipeptide.*?°?°? Thus, substrate binding increases 
pKa of E21 (Fig. 5). Mutations at E21 as well as 
other residues in the polar cluster of YbgH result in 
loss of the transport activity.‘ A very similar cluster 
of polar residues (including the so-called “ExxER” 
motif in TM1) is also observed in the crystal struc- 
ture of the nitrate-H* symporter, NTR1.1, from Ara- 
bidopsis (PDB ID: 40H3).°* Taken together, 
substrate loading-induced increase of pKaou, may be 
a common event during symporter protonation. 


Antiporters 


In contrast to symporters, under condition of limited 
substrate availability, a PMF-driven MFS antiporter 
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would stay in the protonated C;, conformation, wait- 
ing for the binding of a substrate from the cytosol 
side (unless additional regulatory mechanisms 
exist). In an ideal case of the competition between 
substrate binding and protonation in an antiporter, 
the pK, in values of the proton-titratable sites in the 
Ci, state would be higher than pH;, in the absence 
of a substrate, preventing premature deprotonation; 
and substrate loading would decrease pK, in to below 
pH;, thus promoting deprotonation. In the Cout 
state, the pKa out values of the titratable sites would 
be lower than pHout in the presence of the sub- 
strate, preventing premature protonation. Substrate 
releasing would increase pK, our to above pHout thus 
promoting protonation. In short, in the case of an 
antiporter, substrate binding destabilizes the proton 
binding. 

How, then, is such a mechanistic scheme imple- 
mented in an antiporter structure? In the Cm state, 
the titratable residues may be buried in a hydropho- 
bic or electronegative environment to maintain a 
protonation state in the absence of a substrate. Sub- 
strate loading may make them either more accessi- 
ble to hydrophilic groups or more electropositive; 
either way it would promote deprotonation. For 
example, in the multidrug-resistance antiporter, 
E. coli MdfA, D34™?* is a potential titration site in 
the central cavity.” In the Cm crystal structure of 
the MdfA-chloramphenicol complex, D34 is directly 
involved in substrate binding via H-bonds." In addi- 
tion, it has been suggested that loading of a sub- 
strate may decrease pK,in of this titratable site 
through interactions with the positively charged 
motif-B (see below). In another multidrug-resistance 
antiporter LmrP, the membrane-embedded acidic 
residue E327™1°4 was found to be essential for 
high-affinity binding with the cationic substrate 
Hoechst-33342 and for proton transport.°° Upon 
interaction with substrates, this residue is believed 
to become more accessible to solvent, thus easier to 
get deprotonated.*” In E. coli YajR, the H225™7*_ 
E3207 pair in the C-domain is the only titrata- 
ble site inside the central cavity. In the Cout state, 
this pair is away from the basic residue R108™4? 
(i.e., the motif-B in the N-domain).“* In the putative 
Cm state, the H225—-E320 pair is likely to approach 
closer to R108, presumably allowing substrate to 
trigger deprotonation. In general, because of differ- 
ent pK,in values of the titratable sites inside the 
central cavity, substrates carrying different charges 
may induce differential deprotonation at distinct 
titratable sites. This might explain the observation 
that the efflux of different substrates is coupled to 
an influx of different stoichiometry of protons, with 
a net inward movement of +1 electric charge in 
each electrogenic transport cycle.°&°* The membrane 
potential A¥ can only be a driving force for sub- 
strate efflux if each excreted cationic substrate has a 


1568 = PROTEINSCIENCE.ORG 


lower charge than the number of protons translo- 
cated.” Such antiporters require multiple protona- 
tion sites that have pK,in higher than pHm in the 
absence of substrates, but can lower the pK, in val- 
ues in response to substrate binding. Taken 
together, a substrate loading-induced decrease of 
PKa,in may be a common event in the transport cycle 
of PMF-driven antiporters. 


Uniporters and electroneutral transport 
Uniporters mediate facilitated diffusion and, by defi- 
nition, do not require proton translocation across the 
membrane as their driving force. Nevertheless, they 
may also utilize the negative-inside membrane 
potential if either the substrate influx carries posi- 
tive charges or the substrate efflux carries negative 
charges. Although some MFS transporters are 
reported to carry out only uniport processes (e.g., 
GLUT1-4), both PMF-driven symporters and anti- 
porters may behave like uniporters if under a strong 
concentration gradient of the substrate. For exam- 
ple, LacY was shown to carry out both symport and 
facilitated diffusion.” Furthermore, many cationic 
compounds can be either exported or imported 
through the same MFS “antiporter,” depending on 
the direction and strength of the concentration gra- 
dient of the substrate between the two sides of the 
membrane.” For example, LmrP may carry out both 
antiport and facilitated diffusion of cationic sub- 
strates.°? Also, Candida albicans Mdrlp usually 
expels a broad spectrum of electroneutral or cationic 
drugs in a PMF-driven manner; however, it has 
recently been shown that Mdr1p also uptakes a class 
of positively charged isoquinoline derivatives.°° Such 
uptake may be driven by AY via facilitated diffu- 
sion. Therefore, there appears to be no unique struc- 
tural features that distinguish uniporters from 
either symporters or antiporters, except that in uni- 
porters, substrate binding or release is not coupled 
in an obligatory manner to a particular protonation 
status. 

Not all MFS transporters are driven by mem- 
brane potential-mediated electrostatic interaction. If 
the net charge change in the transport cycle is zero, 
the transport is called electroneutral [Fig. 3(C)]. An 
electroneutral transport is often defined by the fact 
as to whether the given transport process is inde- 
pendent of AY. In such cases, another favorable 
chemical potential has to serve as the thermody- 
namic driving force. However, since the presence of 
a membrane potential is a universal feature of all 
living cells, electrostatic interactions may play cer- 
tain roles in affecting the kinetics of the transport 
process, even for electroneutral transport. For exam- 
ple, it would require at least a 55-fold concentration 
gradient (RTIn55= 4 RT) for a monovalent charged 
substrate to move against a 100-mV membrane 
potential, and such a gradient would be unusual 
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Figure 6. Active conformation of motif A of E. coli YajR. A: Side chains in motif A and Asp126 from the charge-relay triad are 
shown as stick models. Ca atoms of conserved Gly residues are shown as spheres. TM11 is colored from cyan at the N- 
terminal end to pink at the C-terminal end. B: Putative role of the inter-domain linker in regulating motif-A. 


under physiological conditions. Therefore, the free- 
energy required for moving a charged substrate 
against the membrane potential is most likely com- 
pensated by an accompanied favorable charge move- 
ment along the membrane potential of either 
protons or ions. Electroneutral symport (i.e., H*/A™ 
co-transport, where A` stands for a monovalent ani- 
onic substrate) has been briefly discussed before,” 
for example, for uptake of some acidic substrates. 
The fungal H*-phosphate (P;) symporter PiPT is 
one such example (PDB ID: 4J05).°' Furthermore, 
electroneutral antiport (i.e, H*/B* antiport, where 
B* stands for a monovalent cationic substrate) is 
more commonly observed in multidrug exporters 
(e.g., MdfA**). In such cases, the energy barrier 
from the Cy, state back to the Cout state may be low, 
such that in the absence of net charge change, a 
negative AGp(s) is sufficient to stabilize the Cout 
state, thus thermodynamically favoring the Cj,-to- 
Cout conformational change [see Fig. 3(C)]. 

As discussed in the Appendix A [Eq. (A4)], in 
the absence of an effect of a membrane potential 
(e.g., in electroneutral antiport), conformational 
changes between Cm and Cou states would be 
driven by the sum of differential-binding energies of 
the transported substances (i.e., AGp(H*) +AGp(S)). 
These energy terms of differential binding are 
intrinsic properties of the transporter, but are not a 
function of substance concentrations. Their sum pre- 
dicts whether influx or efflux of the substrate is 
thermodynamically favored at the steps of conforma- 
tional changes, in the absence of effects of mem- 
brane potential. Nevertheless, since a membrane 
potential determines the equilibrium positions of 
both Cm and Cout states, AGp(H*) and AGp(S) per 
se may be affected by AY, even in electroneutral 
transport. In addition, since different substrates 
may differ in their differential-binding energies 
between the Cm and Cout states, their transport 
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mechanisms in a given transporter are not necessar- 
ily the same. This is particularly true in the absence 
of a membrane potential, whereupon differential- 
binding energy becomes a major energy source. 


Conserved Motifs 

Despite the drastically different sizes, shapes, and 
chemical properties of their substrates as well as dif- 
ferent protonation mechanisms, MFS transporters 
often contain several conserved motifs in addition to 
the 12-TM helix overall structure.°? Among them, 
motifs-A and -B are the most common. Such con- 
served motifs are likely to provide the structural 
basis of common functions of MFS transporters, for 
example, stabilizing one of the two major conforma- 
tions and regulating or responding to either protona- 
tion or deprotonation. 


Motif-A 


3D motif-A. Motif-A (“GxxxD‘*RxGRR’) is located 
in the intracellular loop L2-3 where it connects TMs 
2 and 3, and is the most conserved motif in MFS 
transporters.° For example, of 77 MFS proteins iden- 
tified from the E. coli genome, 62 contain a putative 
motif-A (Supporting Information Table S4 in Ref. 14). 
Although most available crystal structures of MFS 
transporters are determined in the Cyn state, which is 
a non-functional state for motif-A, mutagenesis anal- 
yses have repeatedly shown that motif-A is essential 
for the transport activity in many MFS transport- 
ers.172142,63-71 More importantly, in the Cout crystal 
structure of YajR/3WDO,™ the motif-A is captured in 
its functional state [Fig. 6(A)]. This structure shows 
that the acidic residue D(+5) of the motif-A partici- 
pates in an inter-domain charge-helix dipole interac- 
tion with the N-terminal end of TM11 from the C- 
domain. Other physical components of the motif-A 
and associated structural features include the follow- 
ing: (i) G(+1) interacts with conserved Gly residues 
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from the N-terminal end of TM11 to ensure a tight 
packing between TMs 2 and 11 in the Cou state [Fig. 
6(A)]. (ii) D(+5), R(+9), and a conserved acidic resi- 
due from the C-terminal end of TM4 (i.e., TM4.5 
region) form a charge-relay triad. (iii) G(+8) is 
required for an S-shaped turn of the L2-3 loop. (iv) 
R(+10) serves as the C-cap for TM6. (v) A conserved 
“PET(S)” motif at the loop C-terminal to TM6 stabil- 
izes both R(+10) (using the mainchain carbonyl oxy- 
gen of the Pro residue) and the TM4.5 acidic residue 
of the charge-relay triad (using the sidechain 
hydroxyl group of Thr/Ser), and serves as an N-cap 
for TM8 (using the sidechain carboxylate group of the 
Glu residue) (see the crystal structures of MdfA/ 
4ZOW, YajR/3WDO, and XylE/4GC0). (vi) A positively 
charged residue (prior to the TM11 N-terminal Gly) 
forms a salt-bridge with the TM4.5 acidic residue of 
the charge-relay triad in the Cout state. Most of these 
residues have been shown in mutagenesis analyses to 
play important roles in substrate transport. For 
example, a mutation at Asp128'™*° significantly 
reduces the transport activity of LmrP.®* Mutations 
in the E. coli antiporter MdfA (e.g., S133F and 
A191T/V) that potentially disrupt the PET(S) motif 
affect the transport activity.” Thus, the motif-A and 
its surrounding residues form a conserved “3D motif- 
A.” One important role played by this 3D motif-A is 
to regulate stability of the Cout state. 


Amphipathic helix. Motif-A stabilizes the Cout 
state. To initiate the Co,-to-Cy, conformational 
change, the Cout state must be destabilized, which 
occurs in response to protonation inside the central 
cavity. How then does this 3D motif-A on the cyto- 
solic side of the transporter sense and respond to 
protonation inside the central cavity? In nearly all 
known crystal structures of motif-A containing, 
PMF-driven, MFS transporters, the charge-relay 
triad is always accompanied by an inter-domain pep- 
tide (i.e., L6-7). These linker peptides contain a 
number of polar residues (often positively charged) 
as well as an amphipathic a-helix (a6—7) (Support- 
ing Information Fig. S2 in Ref. 14), the latter of 
which is likely to bind to the inner leaflet of the 
membrane. If this amphipathic a-helix maintains 
membrane association during an inward shift of the 
TM core, the inter-domain peptide will be stretched, 
resulting in a linker movement that likely interrupts 
the charge-relay triad of the motif-A [Fig. 6(B)]. For 
example, in the Cm crystal structure of NRT1.1/ 
40HB, a lysine residue from the linker peptide is 
found to insert into the charge-relay triad.** In the 
Cin crystal structure of MdfA/4ZOW, D(+5) of motif- 
A is surrounded by positive residues from the motif- 
A as well as L6-7 loop. Moreover, in the crystal 
structures of several bacterial peptide transporters 
(POTs), an inter-domain, TM-helix pair is inserted 
between the two domains of the TM core, replacing 
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a6—7, and is likely to move independently from the 
rest of the TM core.2”*%2 Thus, this helix-pair 
insert may function in a way similar to the amphi- 
pathic helix a6—7 in other 12-TM MFS structures. 
In support of such possibilities, mutations in the 
L6-7 that change the charge distribution or disrupt 
the amphipathic property of the a6-7 result in loss 
of transport activity. “>> Interestingly, in the 
multidrug-resistance antiporter MdfA, mutation of 
S204 at the N-cap of «6-7 to a Pro residue enhanced 
the transport activity.” Therefore, the L6-7 and 
associated amphipathic a6-7 are likely to play the 
roles of sensing protonation and regulating the sta- 
bility of the 3D motif-A. 


Motif-A hypothesis. Together, the above data sup- 
port the “motif-A hypothesis”’*: Protonation-induced, 
inward movement of the TM core relatively to the 
membrane is sensed by the amphipathic helix, «6-7. 
This helix further induces a conformational change 
in the L6-7 loop. Through the charge-relay triad, 
the positive charges of the loop weakens the inter- 
domain, charge-dipole interaction between D(+5) of 
the motif-A and TM11, thus destabilizing the Cout 
state of the MFS transporter. 


A-like motifs. Because of the internal symmetry 
of MFS transporters, motif-A could exist at four 
locations, namely L2-3, L5-6, L8-9, and L11-12. 
They are referred to as A-like motifs (or motifs- 
AM3, -A58 Ale? and -A "?), Among the four 
potential A-like motifs, A? and A18 are located 
on the cytosol side, and A15® and AH”? are located 
on the periplasmic side. Sequence analysis of E. coli 
MFS transporters has confirmed the existence of A- 
like motifs in all four locations.!* In some cases, 
more than one such motif exists in a given MFS 
transporter. For examples, E. coli YbgH contains the 
motifs-A and -A842 Mammalian multidrug- 
resistance antiporter Octl contains two complete 
sets of motif-A at both L2-3 and L8-9 on the cyto- 
plasm side, and so does the fungal H*-phosphate 
symporter, PiPT (PDB ID: 4J05).°! Signature motifs 
of sugar porters [Supporting Information Figs. S1(B) 
and S2 in Ref. 74] are in fact components of the 3D 
motif-A? and -A™**. These two motifs are found 
even in a non-PMF driven nitrate/nitrite exchanger, 
NarK” (note that its homolog NarU is proposed to 
be a (Na*/K*)-(NO;/NO,) symporter’®), indicating 
that the motif-A is not directly involved in proton 
translocation. Mutations in these A-like motifs were 
shown to reduce transport activities.t77>"” Their 
main functions are likely to include stabilizing of 
the Cout or Cy, states. The precise mechanisms that 
regulate these A-like motifs in L5—6, L8—9, and L11- 
12 remain to be elucidated. 
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Figure 7. Structure and mechanism of the motif-B. A: Motif-B in the crystal structure of MdfA/4ZOW. Residues of the motif-B 
and surrounding conserved residues are shown in stick models. Backbone of the N- and C-domain is shown in tubes indicated 
in wheat and red, respectively. B: Schematic diagram of motif-B in regulating deprotonation. 


Motif-B 

Motif-B, “RxxqG” in TM4, is present in many MFS 
transporters, including a number of multidrug- 
resistance antiporters.” According to the PFAM 
database,’ motif-B is the only high frequency site 
that contains a basic residue embedded in the TM 
region among MFS transporters (Supporting Infor- 
mation Fig. S2 in Ref. 17). Of the 77 MFS proteins 
from E. coli, 13 transporters possessing the “RxxQG” 
motif are identified. Nine of them, namely EmrB, 
EmryY, HsrA, MdfA, MdtD, MdtL, MdtM, YajR, and 
YebQ, are annotated to be PMF-driven antiporters, 
and the remaining four (KgtP, ShiA, YdfJ, and 
YhjE) are symporters. This classification indicates 
that the simple presence of a motif-B is neither a 
necessary nor a sufficient condition for an MFS 
transporter to be an antiporter. In addition, in FE. 
coli there are six “R™*?xx[S/T/C/H]G”-containing 
transporters, namely AraJ, Bre, DgoT, EntS, Fsr, 
and MhpT; and there are 14 other “R™*?xxxG”-con- 
taining transporters, including ExuT, GalP, MdtG, 
NanT, NepL, UhpC, XylE, YaaU, YcaD, YhhS, YjhB, 
YjiJ, YjjL, and YnfM. Hypothetically, a motif-B 
might appear in the four cavity-helices (i.e., TMs 1, 
4, 7, and 10), which are related by the internal 
pseudo-symmetry of the MFS transporter. However, 
sequence showed that the “RxxQG” 
sequence motif occurs only in the TM4.2 region. Fur- 
ther, two of the nine “RxxQG”-containing E. coli 
MFS antiporters (EmrB and EmrY) do not contain a 
motif-A. Thus, potential functions of the motif-B 
appear to be independent of the motif-A. More 
importantly, mutations in the motif-B in several 
MFS transporters have been shown to abolish all 
transport activity.1%°’1795! Together, these data 
suggest that motif-B plays a common functional role 
in motif-B containing MFS transporters. 


analysis 


MdfA and YajR crystal structures. Among cur- 
rently available crystal structures of MFS transport- 
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ers, the canonical motif-B is only observed in Æ. coli 
YajR and MdfA.1” Similar to motif-A, motif-B is also 
accompanied by a number of residues that are con- 
served in the 3D structure but not contiguous in the 
primary sequence [Fig. 7(A)]. Together, they form a 
“3D motif-B” that contains an extended H-bond 
network. 

Motif-B is unlikely to be directly involved in the 
inter-domain conformational change, because the 3D 
motif-B is embedded entirely in the N-domain. 
G(+5)™+3 residue of motif-B is more conserved 
than R(+1)™*?, even in MFS transporters carrying 
no motif-B, and allows a tight packing of TM4 with 
TM2. Thus, R(+1) appears to be the signature resi- 
due of motif-B. It carries a buried positive charge 
and breaks the H-bonding pattern of the neighbor- 
ing TM1 helix.!” Both features are energetically 
costly, implicating important functions of the “3D 
motif-B.” Mutagenic analysis on MdfA shows that, 
among point mutations at R(+1), only the R(+1)H 
variant confers resistance to most antibiotics tested, 
the R(+1)K variant partial drug- 
resistance activities, whereas R(+1)M/Q/E variants 
lose all activities.°° Results of drug-resistance assays 
on more extensive mutations in the 3D motif-B are 
in agreement with its importance in the transport 
activity." 

In the Cout state, levels 1 and 2 of the cavity- 
and rocker-helices are accessible to solvents from 
the periplasmic side, thus the electrostatic field radi- 
ating from the motif-B is attenuated by the solvent, 
which is of a high dielectric constant [Fig. 7(B)]. In 
the Cm state, in contrast, levels 1 and 2 of the TM 
helices become buried inside the inter-domain inter- 
face, thus the positive electrostatic field of the motif- 
B is strengthened relatively to the Cout state. In 
particular, the level 2 region would become more 
electropositive, thus potentially being able to pro- 
mote deprotonation in the central cavity. For exam- 
ple, in the Cm crystal structure of MdfA/4ZOW [Fig. 
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7(A)], R112*” of the motif-B is 8 A away from the 
protonation site D34™1? in the N-domain. In com- 
parison, in the putative Cj, state of YajR, the titrat- 
able H225™73_—320™!03 pair in the C-domain 
becomes closer to the motif-B than in the Cout state 
of YajR/3WDO, thus favoring deprotonation in the 
Ci, state. Therefore, effects of motif-B do not depend 
strictly on the location of the proton-titratable sites 
inside the central cavity. Furthermore, in 
multidrug-resistance antiporters, many of which 
contain motif-B, loading of substrates (which are 
often small hydrophobic molecules) from the cyto- 
solic side may further strengthen the electrostatic 
field, namely by reducing the dielectric constant 
inside the cavity and/or by bringing in a positive 
charge into the cavity, thus promoting deprotona- 
tion.!” Taken together, in the Cm state the positive 
electrostatic field radiating from motif-B seems to be 
essential for regulating the protonation status of 
motif-B containing transporters. 


B-like motifs. In members of the sugar porter 
(SP) subfamily of MFS, motif-B changes to “RxxxG” 
(where x stands for a hydrophobic residue).’* Some 
of SP transporters (e.g., GLUT1-4) are known as 
uniporters, while others (e.g., XylE) are PMF-driven 
symporters. In the crystal structure of XylE/4GCO, 
the side chain of Gln residue of the canonical motif- 
B is replaced by two ordered water molecules. The 
XylE structure also shows a pair of acidic residues, 
D27™13_F206™® in the vicinity of R133™*? of 
motif-B (“R'"IIGG”), whereas those uniporters, 
instead, have electroneutral polar residues replacing 
the acidic pair. Therefore, the acidic pair in XylE is 
likely to be the potential protonation site, and its 
protonation status is probably influenced by R133. 
These acidic residues have recently been shown to 
be essential for the transport activity of XylE.5? In 
addition, this acidic-residue pair is located closely to 
the central cavity. By expelling solvent, substrate 
binding in the Cout state increases the electronega- 
tivity of the acidic-residue pair, promoting protona- 
tion of the latter. Such a mechanism is another 
example of coupling protonation with substrate bind- 
ing in symporters. Once in the Cm state, R133 
decreases the electronegativity of the micro- 
environment of the acidic residue pair, promoting its 
deprotonation. 

In the fungal H*-phosphate (P;) symporter 
PiPT, motif-B assumes a form of “R!°°VFLG”. In its 
Cm crystal structure (PDB ID: 4J05),°" the basic res- 
idue R139™*? is stabilized by Q1197™?. Two acidic 
residues, D45™!:3 and D48™" are located close to 
motif-B, and this pair is likely to be the protonation 
site. In the Cou state, loading of the negative 
charged substrate into the cavity increases the elec- 
tronegativity of the acidic-residue pair, promoting 
protonation. In the Cm state, the combination of 
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enhanced positive electrostatic field of the motif-B 
and release of the negatively charged substrate pro- 
motes deprotonation. 

Although the sugar-H* symporter, E. coli LacY 
(PDB ID: 1PV7) does not possess a motif-B, a basic 
residue, R302™»?, is buried in the C-domain and 
located at the membrane-level similar to R(+1) of 
motif-B in the N-domain (e.g., R112™*? of MdfA). 
This R302 is located near the potential protonation 
site, H3227™°.2_—325™103 in LacY. Mutations at 
R302 abolish the proton translocation in LacY,**® and 
this observation may be explained using the mecha- 
nistic insights gained from studies of motif-B. 

In all cases discussed above, the protonation 
site inside the central cavity is located on the cyto- 
solic side of motif-B. Increasing the strength of a 
positive electrostatic field inside the central cavity 
in the Cm state relatively to the Cout state seems to 
be a common mechanism of promoting deprotonation 
in MFS transporters. 


Stability Balance Between the Two 
Conformations 

The rocker-switch mechanism of MFS transporters 
suggests that a balance between the Cou and Cin 
states is essential for the proper function of the 
transporters (see Appendix B). This notion has been 
supported by a number of studies. For example, 
point mutations in E. coli LacY that rescue a defec- 
tive mutation at the D(+5) position of motif-A? on 
the cytosol side are found to be clustered around the 
motif-A"'- on the periplasmic side,” implicating 
that balanced stabilities of the Cou, and Cm states 
are important for the transport function. Moreover, 
the C154™5?G point mutation is an extensively 
studied variant of LacY. This mutation is 
transport-inactive, while maintaining substrate 
binding (with affinity even higher than the WT).”° 
The crystal structure of detergent-solubilized 
C154G-LacY (PDB ID: 1PV7) was among the first 
reported MFS crystal structures’ and was found to 
assume the Cm conformation. It was speculated that 
the C154G mutation stabilized the Cm state; never- 
theless, the WT-LacY crystal structure (PDB ID: 
2V8N) was later reported to be also in the Cy,-state.°4 
In fact, a crystallized conformation only suggests that 
it is a dominant population under the crystallization 
condition (including crystal packing). Moreover, a 
DEER study showed that in solution, C154G was pre- 
dominantly in the Cou state (Figs. 2.3 and 2.5 in 
Ref. 20), while the WT assumed more equal distribu- 
tion between the Cm and Cout states (Figs. 2.2 and 
2.4 in the same Reference). Since substrate binding 
per se is shown to stabilize the Cout state,”° the 
tighter substrate binding in C154G seems to correlate 
with a more predominant Cout conformation in solu- 
tion than WT. Therefore, the lost transport activity of 
C154G is likely to be caused by the reduced stability 
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in the Cj, state. This corollary is in agreement with 
the close location of the mutation site to the inter- 
domain interface in the Cm crystal structure (1PV7). 
Intriguingly, the loss of transport activity of 
periplasmic-side mutation C154G can be rescued by 
secondary mutations on the cytosol side of the trans- 
porter. In particular, rescuing mutations include 
E130™45C, K1317™**C, R134C, and G141™>5c, 
These mutations are potentially involved in the inter- 
domain interface in the Cout state, presumably desta- 
bilizing the latter. Again, balanced stabilities of the 
two states appear to be important for transport activ- 
ity. A similar balance was also shown to be important 
for E. coli POT transporter YbgH.*” Mutation D70A 
at the critical acidic residue of motif-A‘?* on the 
cytosol side results in loss of transport activity, pre- 
sumably by reducing the stability of the Cout state. 
Such an activity loss can be reversed by another loss- 
of-function mutation on the periplasmic side, E163A/ 
R297A, which reduces the Cm state stability by dis- 
rupting the motif-A® and a salt-bridge bond 
observed in the Cm crystal structure of YbgH/4Q65. 
In all three of the above studies, loss-of-function 
mutations on the cytosol side and on the periplasmic 
side can rescue each other. These results strongly 
suggest that the balance in stability between the 
Cout and Cy, states plays a key role in the transport 
activity of MFS transporters. 

It is most likely that adjusting stabilities of the 
two major conformations is a common mechanism to 
regulate the activity of most alternating-access 
transporters, including MFS and non-MFS trans- 
porters. For example, transporters of the anti-acid 
system, for example, Arg*/Agm antiporter of E. coli 
AdiC (a member of the LeuT transporter family), 
preferentially perform the transport at low pH 
rather than neutral pH. According to its crystal 
structure (PDB ID: 3NCY*), a sensor for acidic con- 
dition is likely to be a cytosolic, inter-domain, salt- 
bridge bond, R78-D307, which is accompanied by 
Y74.°° This salt-bridge bond likely stabilizes the 
Cout State, even in the presence of a bound substrate 
Arg” ready to be transported. The interaction of 
Arg’ with membrane potential likely provides the 
driving force to the Cout-to-Cyj, conformational 
change. In contrast, Y74 uses its positively charged 
edge of the sidechain aromatic ring as well as its 
hydroxyl group to decrease the pK, of D307, so that 
the strength of the salt-bridge bond is enhanced at 
neutral pH. However, at low pH D307 becomes pro- 
tonated, and the salt-bridge bond dissolves more eas- 
ily, thus shifting the Co y-Cy, equilibrium toward 
the Cm state. In addition, many MFS transporters 
contain rather large exo-membrane domains (e.g., E. 
coli XylE™ and YajR®’), and such regulatory 
domains likely exert their functions through manip- 
ulating stability of certain states. 
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Concluding Remarks 

In the past decade, structural studies on transport- 
ers have provided significant insights into the struc- 
tural basis of the alternating access mechanism. In 
PMF-driven electrogenic MFS transporters, interac- 
tion of protonation and the negative-inside mem- 
brane potential provides a significant part of the 
driving energy for the Cout-to-Cın conformational 
change. For symporters, substrate loading triggers 
the protonation in the Cout state; while for antiport- 
ers, substrate loading triggers the deprotonation in 
the Cm state. Conserved motifs A and B of MFS 
transporters play important roles in common func- 
tions of MFS. Balanced stabilities between the Cout 
and Cm states are essential for the transport func- 
tion. Future work will be required to define and ver- 
ify details of related mechanisms. 
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APPENDIX A: THERMODYNAMIC DESCRIPTION 
OF MFS TRANSPORTERS 


In a transport process, the “external” Gibbs free energy 
includes terms of Apuy, Auipy+), and Ajs.” The following 
is true for a cycle of PMF-driven transport of (electro- 
genic) substrates: 


Apy + Appr) + Ang =—-@ (A1) 


(< 0. Q is the heat released in one transport cycle); 
where 


def 


Apy = FAY (< 0); (A1.1) 
Aum+ = —2.3RTApH (A1.2) 

(< 0. ApH = pH;,—pHou ~ +0.6); 
Aus, & RTIn((Slp/[Sh,) (41.3) 


(> 0, if the transport is against the substrate gradi- 
ent. L and R stand for loading and releasing states, 
respectively). 

Each term of external energy can be further 
divided into sub-terms (see Fig. 3). 


Ang) =AGy(S) +AGp(S)+AGR(S); (A2) 
AG (S) © -RTIn([S],/KaL) (A2.1) 

(free energy of substrate loading); 
AGp(S) & RTIn(Kar/Kat) (A2.2) 


(substrate-associated intrinsic differential free 
energy of the conformational change from the L to R 
state); 
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def 


AGR(S) £ RTIn((S]p/Kan) (A2.3) 


(free energy of substrate releasing). 
Similarly, the following equations for proton- 
related free energies apply: 


Augrj=AGL(H*) +AGp(H*)+AGR(H*) — ((A8)) 


AG, (Ht) & -2.3RT (pKa out (H+) —pHo.) (A3.1) 


(free energy of protonation in Cout); 


AGp(H*)  2.3RT (pKaout(H*) —pKain(H*)) 


(A3.2) 


(protonation-associated intrinsic differential free 
energy of the Co,;-to-Cy, step); 


AGr(H*) © 2.3RT(pKain(H*)—pHj,) (43.3) 
(free energy of deprotonation in Crn). 
Further, 

AG}., + AGP, 9 =AGp (H*) +AGp(S) (A4) 


(where AG®,_, and AG?_, are free energy changes asso- 
ciated with conformational changes in the absence of 
membrane potential). 


APPENDIX B: KINETIC DESCRIPTION 


A kinetic discussion can be made based on the 
“tightly coupled chemical model.” Based on the 
assumption that, in a symporter, substrate binding 
is coupled with protonation, the following equations 
apply: 


Cout 
( Sluha (Haul? 
out” 1 out’ a) 
OF OR 
kı ka 
ke || k. kz expl- EEH) || kz expl- C2550) 
[H*]ink?-s [S]ink?.4 
() = Gr) = Gs) 
7 ke “(ky 
H+ s 
C, 


n 


AGL(S)=-RTIn((k?[S]) /k-1) (B1) 


(where k? and k-ı are forward and reversal kinetic 
coefficients for substrate loading); 
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AGL (H*)=—RTIn((A2[H*]o,.)/k-2)  (B2) 


(where k9 and k_» are forward and reversal kinetic 
coefficients for proton loading); 
AGQ.1=—-RTIn(k3/k-3) (B3) 


(where kg and k-3 are forward and reversal kinetic 


coefficients for the Coutto-Cm conformational 
change in the absence of Apy); 
AGr(S)=-RTIn(k4/ (k? 4[S]R)) (B4) 


(where k4 and k°, are forward and reversal kinetic 
coefficients for substrate releasing.); 


AGR(H*)=—RTIn (ks /(k° 5 [H] n) (B5) 


(where ks and k°; are forward and reversal kinetic 
coefficients for proton releasing.); 


AG? 9 =—RTIn(ke/k-6) (B6) 


(where kg and k-e are forward and reversal kinetic 
coefficients for the Cm-to-Cout conformational change). 
In addition, 


AGp (H*)=RTIn((k9ks) / (R-2k° 5) (B8) 
and 
(ki kgkskakske) /(k-1k-2k-3ho yk? sk-6) = 1 (B9) 


(i.e., the tightly coupled chemical model; equivalent 
to Eqs. (A4)). 

Similarly, based on the assumption that in an anti- 
porter substrate binding is coupled with deprotonation, 
one has the following diagram as well as equations simi- 
lar to the above discussion for the symporter. 


Cout 
Co owas =) 
“TSloutla 
ke || k.s k.s expl- EY) || kz exp C=2EAV) 
ni. 
F n = 0) 
Cin 


APPENDIX C: KINETIC OF IN VITRO SUBSTRATE 
BINDING ASSAYS 


Substrate binding is often studied in a solubilized form 
of a transporter. The results from such a binding assay 
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usually represent transporters at different states with 
or without substrate bound. Major differences in 
kinetics of a membrane embedded transporter compared 
to the solubilized form are the lack of external driving 
force and identical substrate concentration for both sub- 
strate loading and releasing in the solubilized form. In 
the following, we assume that the transporter maintains 
the two major conformations in the soluble form as it 
does in the membrane. In addition, the affinity of sub- 
strate binding in each conformation is assumed to be the 
same in both the membrane embedded and the solubi- 
lized form. (This assumption may need to be verified in 
each case.) Each conformation is of distinct properties, 
in particular binding affinity (Ka) of the substrate. Note 
that Kg is related to binding enthalpy, AH, by the van’t 
Hoff equation. The kinetic relationship between these 
states is shown in the following scheme. 


Cout 


substrate free 
ro 
ee 
> 
> 
= 
—E S 
jo 
substrate bound 


Cin 


Equilibration Conditions 


[CoutS]/[Cout] =k1[S]/k-1 (C1.1) 
[CmS]/[CouS] =k2/k-2 (C1.2) 
[Cm]/[CmS] =k3/k? [S] (C1.3) 
[Cout]/[Cm] =k4/k-4 (C1.4) 
(k?koksk4)/(k-ık-2k?;k-4) = (C1.5) 


Although there are seven independent kinetic 
variables in Eq. (C1.5), there are only three inde- 
pendent thermodynamic variables (e.g., equilibrium 
constants). 


Dissociation Constants 


def e_1/k°= [S][Cout] /[CoutS] (C2.1) 


Kas © ks/k®3= [S]|Cm]/[CrmS] (C2.2) 


Kaspp. © [S] (Cou) +{Cul)/((CouS] +[CmS]) (C2.3) 


= Kai/(1 +k2/k-2) + Kag/(1 + k-2/k2)  (C2.4) 
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= Kai(1 +k-4/ka)/(1 +k2/k-2) 
E Ka3(1 +k4/k-4)/(1 +k_2/k) 


(C2.5) 
(C2.6) 


Whilst dissociation constants Kg, and Kas; in 
the Cou and Cy, states, respectively, and their 
related free energy terms (AG,(S) and AGR(S)) are 
of interests, what is directly measurable in in vitro 
binding assays such as SPR, ITC, and SPA is usu- 
ally the apparent dissociation constant of the mix- 
ture system, Kg app. 


Partition Between the Cou and Cm 
Conformations 

Here, we define the ratio of the concentration of 
transporters in both substrate-free and substrate- 
bound Cout states to that of Cm states as conforma- 
tional ratio f, which is a function of [S]. 


f(S]) = (Cou) +[CoutS])/((Cm8] +Cm])  (C3.1) 
= (k-2/k2) (1 +Kas/[S])/(1 +Kas/[S])  (C3.2) 
= (ka/k-4) (1 +[S]/Ka1)/(1 +[S]/Kaz)  (C3.3) 

floc) = k-2/ke (C3.4) 
f(0) = ka/k-4 (C3.5) 
= 1 +K, 
E fear) tan 29 
= f(0) (1 +(F(œ)/F(0)) [S]/Kas)/(1 +18]/Kas) 
(C3.7) 


The coefficients K0) and ftoo) are the Cou-Cin 
conformational ratios in the absence of substrate 
and in super-high concentration of substrate, respec- 
tively. They are important properties of the trans- 
porter for a given substrate, and can be chosen as 
independent thermodynamic variables. Experimen- 
tally, the curve of f([S]), including f(0) and foo), can 
be measured with assays such as smFRET or DEER. 
The third independent variable, for example, Ka; 
can be determined by fitting the observed f[S]) 
curve with the calculated one using Ka, as a vari- 
able [in Eq. (C3.7))]. 

In general, the Gibbs free-energy of conforma- 
tional change from the Cout to Cy, state at the sub- 
strate concentration [S] is given by: 

AG((S]) $ AH-TAS = RTIn(f((S})) (C3.8) 
where the enthalpy AH and entropy AS terms are 
commonly considered to be independent of tempera- 
ture. Further, lets assume, for a given [S], AU[S]) is 
measured at two temperatures, Tı and Tə, with that 
AT T.—T, (> 0), T (T,+T2/2, and AT < T. 
According to the van’t Hoff equation, the following is 
applied: 
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AH = RT\n(fi/f2) (T/AT) (C3.9) 
where fı and fọ are the conformation ratios meas- 
ured at Tı and To, respectively. Therefore, by meas- 
uring f([S]) at two temperatures, one may determine 
AG, AH, and AS at any specified substrate concen- 
tration [S]. 


Distribution of the Transporters in the Four States 


[all] = [Cout] +{CoutS] +CmS] +[Cin] (C4. 
Rx((S]) © [Cout]/[all] (C4.2) 
= 1/((1+1/f((S])) (1+[S]/Ka,)) 

R,([S]) & [CoutS]/fall] (C4.3) 
1/((1+ 1/f([S])) (1+Ka.1/[S})) 
Ra((S]) & [CiaS] /[all] (C4.4) 
=1/((1+f((S])) (1+Kaz/[S])) 
Ra4((S]) Œ [Crm] /[all] (C4.5) 


=1/((1+f((S])) (1+ [S]/Kas)) 


Of Rx(0), Ro(co), R3(co), and R4(0), three of them 
can be chosen as independent variables. The fourth 
one can be determined by the equation Ry+R2+R3+ 
R= 1 [ie., Eq. (C4.1)]. 


Relationship Between Dissociation Constants 
and f([S]) 


Kai=([S) (/F((S))—1/F(co))/(/F(0)-L/A(IS})) (C5.1) 
Kas=(S] (f([S])-F(ce))/(F(0) -£ (IS) (C5.2) 
Kaapp. = Kai/(1 +1/f(co)) + Kas/(1 +f(œ%)) (C5.3) 
= Kaa(l +1/f(0))/(1 +1/F(0o)) (C5.4) 

= Kaa(1 +f(0))/(1 +F(co)) (C5.5) 


It can be proved that Ka, is the substrate con- 
centration at which 1/f([S]) equals (1/f(0) +1/floo))/2, 
that is, the mean of 1/A[S]). Kas is the substrate 
concentration at which f([S]) equals (A0) +ftoo))/2, 
that is, the mean of f([S]). Similarly, Kaapp. is the 
substrate concentration at which f([S]) equals (R,(0) 
+Ro(oo))/(R3(oo) +R4(0)), that is, the ratio of the 
mean fraction of Cout to that of Cin. 


Further, 
Kas/Kai= k{ks/(k-1k°3) = k-2k-4/(kək4) (C5.6) 
=f (cc) /f(0) 
f(0)=(1-Kaapp./Kas)/(Kaapp./Kai-1)  (C5.7) 
f(œ) = (1—Ka3/Kaapp.)/(Ka1/Kaapp.— 1) (C5.8) 
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Differential Free Energy of Substrate Binding 


def 


AGp(S) = RT1n(Ka3/Ka,) (C6.1) 


= RTin(f(co)/f(0)) (C6.2) 

AGp(S), which is determined solely by the ratio 
of foo) to K0), is an important property of the trans- 
porter (a symporter in this case) for a given sub- 
strate. It indicates the thermodynamic preference of 
the substrate movement in the absence of external 
driving forces (e.g., membrane potential). When 
AGp(S)<0, the movement is inward; and when 
AGp(S) > 0, the movement is outward. 


APPENDIX D: HOW TO ESTIMATE HYDROPHOBIC 
MISMATCH FORCES 


The atomic solvation parameter (o) of a membrane 
protein (MP) is estimated to be 


(D1) 
o = 20-25 cal/mol/A” ~ 0.1 kJ/mol/A” 


(see Ref. 89). In the absence of a charge—potential 
interaction, the MP is at an equilibrium position (i.e., 
Ax = 0). Thus the hydrophobic mismatch force is 


At a large displacement, 
falAx) larza ~ OC (D2.2) 


where c is the circumference of the MP (typically 
~100 A). So, the larger the size of the MP is, the 
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stronger the hydrophobic mismatch force may 
become. For a small rigid-body displacement, the 
hydrophobic mismatch force can be written as 


falAx) = (1-exp(—Ax/A))oc (D2.3) 


where 4 indicates the rigidity of the membrane. The 
larger / is, the more sticky (dissipative) of the mem- 
brane—MP interface would be. Typically, å is in the 
order of 10 A. 

In the presence of the membrane potential AY 
(or more precisely of a charge—potential interaction), 
the electrostatic force is 


fe= FAY/d (D3) 


where F is the Faraday constant (typically FAY ~ 
10 kJ/mol) and d is the membrane thickness (~30 
A). Assuming a rigid-body displacement, at the equi- 
librium position the two forces balance each other: 


(1-exp(—Ax/1))oc = FAY /d (D4.1) 
Thus, 
Ax/A=—In(1-FAYW/(dac)) ~ 1/30 (D4.2) 


Note that (FAY/d)Ax would be the electrostatic 
energy of the displacement, and (oc)Ax would be the 
hydrophobic mismatch energy assuming a constant 
maximum force during the displacement. Equation 
(D4.2)) suggests that even for a soft membrane (e.g., 
A 10 A), the displacement Ax caused by electro- 
static interaction is likely to be small (unless the 
MP is not a rigid body, e.g. when there is a large 
conformational change). 
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